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Magnetism in 2D networks of the acetylene-bridged transition metal porphyrins M(P)-2(C-C)-2 
(denoted P-TM), and oxo-TM-porphyrins OM(P)-2(C-C)-2 (denoted P-TMO), is studied with the 
density functional theory (DPT) and the self-interaction corrected pseudopotential scheme (pSIC). 

Addition of oxygen lowers magnetism of P-TMO with respect to the corresponding P-TM for most 
of the first-half 3d-row TMs. In contrast, binding O with the second-half 3d-row TMs or Sc increases 
the magnetic moments. Ferrimagnetism is found for the porphyrin networks with the TMs from 
V to Co and also for these cases with oxygen. This is a long-range effect of the delocalized spin- 
polarization, extended even to the acetylene bridges. 

PACS numbers: 68.35.bm, 68.35.Dv, 68.65.Cd, 75.10.Lp, 75.25.Dk 


I. INTRODUCTION 

Porphyrins are popular molecules in living organisms. 
Their ability to bind the transition metals is utilized in 
biological processes significant for the human-blood func¬ 
tions - employing Fe in heme. It is also responsible for 
the photosynthesis in plants - building a complex with 
Mg, namely chlorofileii The oxo-metal-porphyrins, such 
as oxotitanium porphyrin (P-TiO), can be used for solar- 
activated water splitting - what has been theoretically 
proposed to occur due to low-lying ligand-to-metal in¬ 
tramolecular charge-transfer states^ and recently experi¬ 
mentally confirmedi^ 

Porphyrins or phthalocyanines can be assembled at 
metals^i^ and other surfaces^ forming 2D networksii In 
turn, the surfaces or additional ligands change the charge 
state of the transition metal bound to porphyrins with 
respect to the bare 2D network. This is because different 
valence states of a metal, metal-oxygen moiety, or metal- 
Li complex lay closely in the energy4“— Not only the 
single layers of the metal-organic networks exist, also the 
double- and triple-decker layers of metal-phthalocyanines 
have been grown, on both the ferro- and antiferromag¬ 
netic substratesji^ 

For catalytic reactions with the intramolecular charge 
transfer, such as the water splitting, the intermolecu- 
lar connections must be very weakly conducting (non- 
covalent); for example the hydrogen bridges (0-H...0). 
For the spintronic applications, on the other hand, these 
connections should be conducting. For this reason, we 
have chosen the acethylene bridges. Various choices of 
the intermolecular-bridging type were examined experi¬ 
mentally for the energy transfer rateS )^^d^ and a role of 
the frontier molecular-orbitals was theoretically studied 
for the donor-acceptor dimersi^ 

A new branch of material science grows on top of po¬ 
tential utilization of the 2D organic layersi^ in the elec¬ 
tronic devices; as well as the magnetic nanoparticles de¬ 
posited on various substrates?^ Selective incorporation 
of metal atoms into organic templates enables wide func¬ 
tionalization of the layerSfi^ hence many applications are 


plausible. The 4-fold symmetry of porphyrins and ph¬ 
thalocyanines with the transition metals makes their net¬ 
works similar to the Heusler compounds, which are very 
strong magnetsJ^ There are review articles on the spin¬ 
tronic topics in the 2D networks The photovoltaic 
thin films based on the metal-phthalocyanine 2D net¬ 
works - organized in the AA-stacking order - conduct the 
current very well in the columnar directionj^^ Similar ef¬ 
fect should be observed in metal-porphyrins. For tech¬ 
nological reasons, the conductive properties of the single 
metal-porphyrin molecules between the electrodes?^ and 
the ID metal-porphyrin chains^! were also investigated. 
Recently, high-mobility field-effect transistors have been 
constructed by building the ABAB-type stacking of the 
phthalocyanine-VO and -TiO layers?^ 

In this work, we focus on the covalent 2D networks of 
porphyrins with all Sd-row transition metals and TMO, 
bridged with the acetylene moiety. The model struc¬ 
tures of our systems are presented in Fig. 1. We want to 
obtain a material with strong and long-range magnetic- 
exchange interactions. Thus the covalent intermolecular 
connection was chosen, since it has been found that the 
covalent inorganic materials - as hosts for the transition 
metal impurities - are strongly magnetically coupled for 
long distances 


FIG. 1. 2D network of porphyrin-MnO with acetylene 
bridges; top view and side view. 

Interestingly, another molecular system - similar to 




2 


that studied in this work, but assembled at Cu and 
containing Fe - was demonstrated to switch the easy 
magnetization-direction after adding the oxygen.— We 
report the ferrimagnetism, similar to that theoreti¬ 
cally predicted by other authors in porphyrins with Mn 
and connected by 4-bromophenyl to form ID magnetic 
chains This effect can be mapped on the atomic scale 
by the measurement of magnetic resonance spectra at gi¬ 
gahertz frequencies using X-ray magnetic circular dichro- 
ism (XMCD).^ 


II. THEORETICAL DETAILS 

We performed the density-functional theory (DFT) 
calculations using the Quantum ESPRESSO (QE) 
suite of codes,— which employs the plane-wave basis 
set and the pseudopotentials to describe the core elec¬ 
trons. The exchange-correlation functional was cho¬ 
sen for the gradient corrected Perdew-Burke-Ernzerhof 
(PBE) parametrization.— The ultrasoft pseudopotentials 
were used with the energy cutoffs 35 Ry and 400 Ry 
for the plane-waves and the density, respectively. The 
Monkhorst-Pack uniform k-mesh in the Brillouin zone 
(BZ) has been set to 10 x 10 x 1, and the Fermi-surface 
energy broadening parameter 0.02 Ry was chosen for 
a better convergence. The vacuum separation between 
the periodic slabs was 40 A. The pseudopotentials for 
the atoms from Sc to Co were modeled with the semi¬ 
core states (3s and 3p) included in the valence band, in 
contrast, the Cu and Zn valence shells were constructed 
without the deeper states. 

We started with the geometry optimization of a sin¬ 
gle molecule, namely porphyrin-Ti, using the B3LYP 
method which is equivalent to the DFT scheme of the 
BLYP-type with 20% of the exact exchange (EXX) 

This step was done with the quantum chemistry package 
TURBOMOLE,— which represents the wavefuntions in 
the gaussian basis set; we used the correlation-consistent 
valence double-zeta atomic basis set with polarization 
functions for all atoms (cc-pVDZ)— Obtained geome¬ 
try and the lattice vectors, derived from the molecular 
size, were inserted as an input for the geometry optimiza¬ 
tion of the periodic structures, calculated by means of the 
generalized gradient approximation (GGA) scheme (with 
the PBE parametrization) using the QE tool. The above 
procedure is accurate enough, since the central ’squares’ 
of porphyrin and phthalocyanine derivatives do not relax 
much after binding a metal— 

It is known, that the DFT approach underestimates 
the energy gaps and the energetic positions, with respect 
to the Fermi level, of the localized d- or f-shells. This 
fact has consequences in the description of magnetiza¬ 
tion. In order to improve the treatment with the lack of 
the exact exchange, we applied the pseudopotential self¬ 
interaction correction (pSIC) method proposed by Filip- 
petti and Spaldin;^ and implemented by us in the QE 
package— The pSIC method is superior to the DFT-I-U 


approach^ for two reasons: (i) the correction is param¬ 
eter free, unlike the DFT-I-U parameters: the Coulomb 
U and the exchange J, (ii) the correction is applied to 
all atomic shells, not only d- or f-shell of the transition 
metals or the rare earths. The pSIC-kernel includes the 
Hartree and the exchange-correlation potentials, 
calculated on the orbital-density nf{r), dependent on 
spin cr, and built using the atomic pseudopotential or¬ 
bital, The main equations are: 

I[<(!•)]) 

<(!■) =vt 

nk 

The occupation numbers pf are obtained from the pro¬ 
jection of the Kohn-Sham states onto the pseudopo¬ 
tential atomic-orbitals ipi (r) and are the Fermi-Dirac 
occupations. Usefulness of this method has been demon¬ 
strated in a number of various applications— r— For this 
work, the most important are the relations between the 
energetic positions of the shells: 3d-TM, 4s-TM and 4p- 
TM, 2p-0 and 2p-N. Although, we included the SIC for 
all atomic shells, these corrections are ’on-site’ - the same 
like in the DFT-I-U approach and contrary to the Na- 
gaoka model, which includes the ’inter-site’ strong corre¬ 
lations on the parametric grounds . Using the pSIC 
approach instead of the DFT-I-U scheme is especially jus¬ 
tified for the molecular crystals and 2D metal-organic 
frameworks, due to the fact they possess the fiat band 
structures— In constrast, the ordinary strongly corre¬ 
lated systems such as the metal oxides and diluted mag¬ 
netic semiconductors possess only d-type or /-type flat 
bands and the rest of projected atomic states have usual 
semiconducting bandwidths. 


III. RESULTS AND DISCUSSION 

In this study, we focus our attention on magnetism, 
which is a phenomenon sensitive to the system geometry. 
As mentioned in the theoretical details section, the ge¬ 
ometry was optimized for all studied cases. Positions of 
transition metal atoms above the plane of the nitrogens 
square and the oxygen-TM bond lengths are collected 
in the Table SI in the supplemental information- 
inspecting the atomic gradients varied during the geom¬ 
etry optimization and the corresponding magnetizations 
of the systems, we noticed that a tiny change of the 
TM-positions or the TM-0 bond-lengths can influence 
magnetizations. This effect is similar to the phenomenon 
observed in the perovskites^ and also for metal atoms 
at graphene and graphite— In our cases, however, the 
geometric effect does not alter main trends in the results 
reported below. 
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FIG. 2. Total and absolute magnetizations (in ^ b ) for metal-porphyrins and oxo-metal-porphyrins with the acethylene bridges, 
obtained with the GGA and pSIG approaches using the GGA-optimized geometry. 


A. Magnetizations 


In Fig. 2, the total and absolute magnetizations in the 
elementary cells, obtained within the GGA and the pSIC 
frameworks, are presented for porphyrins with all TM 
and TMO additions. The total magnetization, defined 
via the up and down spin-densities, n^(r) and U4,(r), as 


My 


total — 


j n-^{r) dr - J n^(r) dr, 


reflects the summed magnetization of the sample seen 
from a distance. In contrast, the absolute magnetization, 
defined as 

Mabsoiute = J I U^(r) - n|(r) | dr, 


means a sum of the local magnetizations regardless their 
signum. Large difference between the total and the ab¬ 
solute magnetizations gives an information on the space 
separation of the alterred magnetic moments, i.e. the 
ferrimagnetic character of the sample. 

Ferrimagnetism is well pronounced for the porphyrins 
with metals from V to Co, and it is the strongest in 
the Mn, Fe and FeO cases for both calculation meth¬ 
ods, GGA and pSIC, and additionally in the V and VO 
cases obtained with the pSIC scheme. The spin-density 
map for the porphyrin-Mn network is presented in Fig. 3 
for the pSIC method. The vanadium spin-density map 
is similar and included in Fig. S2 in the supplemental 
information.— It is clear that the nitrogen atoms are 
spin-polarized in the TM vicinity. 

The effect of oxygen addition on magnetization is very 
strong. In the case of the TM atoms from the first half 


of the 3d-row, the oxygen addition strongly reduces the 
magnetic moments; except the porphyrin-ScO network. 
In contrast, for the TM atoms from the second half of 
the 3d-row, addition of oxygen causes an increase of the 
magnetic moments. This is an interesting result - similar 
to the results for co-doping of metal-phthalocyanines 
with Liii - and might be utilized in spintronic devices. 



FIG. 3. Spin-density map of porphyrin-Mn, obtained with 
the pSIC approach. 

To get further insight into the electronic structure and 
magnetism, it is useful to examine the Lbwdin popula¬ 
tion numbers. The differences between the Lbwdin pop¬ 
ulation numbers for the spin up and spin down are the 
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TABLE I. Lowdin spin-polarizations (in fis) for the 2D porphyrin networks with TMO (upper part) and TM (lower part); 
calculated with the pSIC scheme. Different carbon atoms are indexed as follows: Cl denotes carbons adjacent to N, C2 
carbons terminated with H, C3 carbons adjacent to acetylene, and C4 carbons of acetylene - they are presented in Fig. 4. The 
corresponding GGA numbers are given in parenthesis. 


shell 

Sc 

Ti 

V 

Gr 

Mn 

Fe 

Go 

Ni 

Gu 

Zn 

O sp 

0.720 

0.0 

-1.570 

-0.215 

0.043 

1.670 

2.265 

1.546 

1.636 

1.491 


(0.577) 

(0.0) 

(-0.149) 

(-0.062) 

(0.577) 

(0.728) 

(1.392) 

(1.403) 

(1.435) 

(1.389) 

TM d 

0.003 

0.0 

0.940 

2.036 

2.991 

2.612 

3.026 

1.562 

0.506 

0.001 


(-0.016) 

(0.0) 

(1.154) 

(1.684) 

(2.424) 

(1.128) 

(1.579) 

(1.749) 

(0.677) 

(0.025) 

TM sp 

0.007 

0.0 

0.005^ 

0.038 

0.087 

0.015 

0.008 

-0.005 

-0.044 

-0.006 


(-0.001) 

(0.0) 

(0.020) 

(0.021) 

(0.024) 

(0.006) 

(0.013) 

(-0.010) 

(-0.043) 

(-0.009) 

N sp 

0.000 

0.0 

-0.270 

-0.047 

-0.036 

0.051 

0.057 

0.101 

0.158 

0.008 


(0.004) 

(0.0) 

(-0.019) 

(-0.032) 

(-0.030) 

(-0.007) 

(0.003) 

(0.135) 

(0.164) 

(0.011) 

G1 sp 

0.000 

0.0 

-0.120 

0.024 

0.006 

0.109 

0.117 

-0.005 

-0.007 

-0.001 

G2 sp 

0.000 

0.0 

-0.057 

0.007 

0.001 

0.056 

0.047 

0.007 

0.007 

0.001 

G3 sp 

0.003 

0.0 

0.052 

0.003 

-0.005 

-0.051 

-0.057 

-0.004 

0.007 

0.001 

G4 sp 

0.001 

0.0 

0.006 

0.003 

-0.001 

-0.006 

-0.003 

-0.003 

0.001 

0.001 

TM d 

0.0 

1.570 

1.630 

3.820 

4.158 

3.940 

3.046 

0.0 

0.420 

0.0 


(0.0) 

(0.989) 

(2.510) 

(3.656) 

(3.543) 

(2.092) 

(-1.044) 

(0.0) 

(0.490) 

(0.0) 

TM sp 

0.0 

0.014 

-0.971’^ 

0.140 

-0.235’* 

0.074 

0.046 

0.0 

-0.027 

0.0 


(0.0) 

(0.199) 

(0.197) 

(0.163) 

(0.146) 

(0.013) 

(-0.092) 

(0.0) 

(-0.018) 

(0.0) 

N sp 

0.0 

-0.063 

-0.290 

-0.074 

-0.470 

-0.109 

0.062 

0.0 

0.136 

0.0 


(0.0) 

(-0.019) 

(-0.057) 

(-0.068) 

(-0.044) 

(-0.036) 

(0.011) 

(0.0) 

(0.129) 

(0.0) 

G1 sp 

0.0 

0.057 

0.210 

0.034 

0.093 

0.172 

0.109 

0.0 

-0.004 

0.0 

G2 sp 

0.0 

0.027 

-0.019 

0.010 

0.009 

0.050 

0.079 

0.0 

0.005 

0.0 

G3 sp 

0.0 

0.002 

-0.002 

-0.009 

-0.039 

-0.071 

-0.056 

0.0 

0.002 

0.0 

G4 sp 

0.0 

0.003 

-0.095 

-0.004 

-0.004 

-0.004 

-0.007 

0.0 

0.000 

0.0 

^ 0.032(s), 

-0.27(p) 



^ 0.006(s), 

-0.977(p) 



0.099(s), 

-0.334(p) 




FIG. 4. The indices for carbon atoms used in the Table 1. 


spin polarizations of chosen atomic orbitals. They are 
listed in Table 1, for all porphyrins with TM and TMO. 
We present the numbers obtained mainly with the pSIC 
scheme. The corresponding GGA results are also listed, 
for a comparison, for the 2p-states of the oxygen and 
nitrogen atoms and the TM atomic orbitals. 

The pSIC method usually tends to localize the 
electrons^ and delocalize the holes.— With the DFT-bU 
method, the 3c?-TM orbital occupations often increase 
and the sp-orbital occupations of the TM-neighbours de¬ 
crease, with respect to the DFT result. For the pSIC 
method this is also usually true, due to the fact that 
the self-interaction correction is stronger for the localized 


shells. Netherveless, this is not a ’golden rule’ and for 
some cases it might be broken. The cases of porphyrin-V 
and -VO are special. Without oxygen, the sp-orbitals of 
V are substantially occupied, because they overlap much 
with the neighbouring nitrogens. It is interesting that 
with the pSIC aproach, the polarization of sp-Y is anti¬ 
ferromagnetic (AF) to Sd-V. Large spin-polarizarions of 
nitrogens - around -0.3 pB, which gives -1.2 ps in total - 
are not surprising. Nitrogen is an element which usually 
couples antiferromagnetically to TM - as for instance in 
(GaN,Mn) dilute magnetic semiconductor— Moreover, 
the sp-Y orbitals overlap and couple ferromagnetically 
(FM) with the sp-N shells. The fact that the 4p-shell 
of V in porphyrins is substantially occupied can be com¬ 
pared to the similar effect observed in (GaAs,Mn)^ and 
Si:Mn— On the other hand, in the GGA approach, the 
sp-Y shell couples ferromagnetically to 3d-V, since the 
strong SI error is not corrected for the sp-N orbitals. 

For the porphyrin-VO, the sp-N is also AF-coupled 
with respect to 3d-V. With the difference that the strong 
coupling of sp-Y with N-sp is replaced by the V inter¬ 
action with oxygen, which is more strongly polarizable 
than nitrogen. The spin polarization of oxygen is around 
-1.6 pb and antiferromagnetic with respect to 3d-Y. The 
sp-N localized spin is -0.27 pB, which is similar to that 
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for porphyrin-V. However, the largest spin polarization of 
the sp-N orbitals is in the case of porphyrin-Mn, around 
-0.47 pb- Interestingly, the coupling of sp-N with re¬ 
spect to 3(i-TM is antiferromagnetic not for all studied 
2D networks - the exceptions are: porphyrin-Co, -Cu, 
-FeO, -CoO, -NiO and -CuO. 

This is a consequence of the 4-fold crystal symmetry 
and the electronic filling of the 3d-TM orbitals. The TM 
atoms possess valence electrons of the s— and d—type 
and need to make the bonds with 4 N atoms and oxygen. 
Formation of the local magnetic moment at the TM atom 
means less electrons involved in the chemical bonds, and 
a hole might be delocalized over these bonds. The lack of 
TM electrons with the spin up at the TM-N bonds - for 
the cases with the d-shell less occupied than it is for Mn 
- causes the domination of the spin down polarization at 
these bonds, and the magnetic coupling of N atoms to the 
TM is antiferromagnetic. The coupling changes signum 
with the growing number of electrons at the d —shell of 
TM, and the cases with Ni and Cu are ferromagnetically 
oriented to the N atoms. Large spin polarization of the 
sp-0 orbitals is very promissing for tailoring spintronic 
devices. It is worth to notice that the orientation of the 
spin coupling at oxygen with respect to 3d-TM is usually 
the same as that of the coupling between sp-N and 3d- 
TM. 

Let us have a look at carbons of the porphyrin 
molecules and bridges; these adjacent to nitrogens, de¬ 
noted Cl, and farther, denoted: C2, C3 and C4 - see the 
caption of Table 1 and Fig. 4. The carbon atom denoted 
Cl bears a spin which couples AF to 3d-V (with local Cl 
moment of -0.12 ps) for the VO case, and FM (0.21 ps) 
for the porphyrin-V. Other spin polarizations at Cl also 
deserve our attention: for porphyrin-Fe (0.17 ps), P-Mn 
(0.09 pb), P-FeO and P-Co (0.11 pb), and P-CoO (0.12 
Pb). The carbon atoms denoted C2 polarize mostly in 
the case of porphyrin-Co (0.08 pb)- Even far from the 
TM atom, C3 atoms, are spin-polarized around 0.04-0.07 
Pb for porphyrin-Mn, -Fe, -FeO, -Co, -CoO. Remarkably, 
in porphyrin-V, the carbon atom of the acetylenic-bridge, 
namely C4, couples antiferromagnetically to Cl, and it 
is spin-polarized of about -0.1 pB. 

The results presented in this subsection align with the 
experimental and theoretical conclusions of the work on 
porphyrin-Cr adsorbed at the Co surfacei^ Namely, the 
local spin-moments induced at nitrogens are AF-coupled 
to Cr, and the spin of less than halfly-occupied 3(i-shell 
of Cr is AF-coupled to the more than halfly-occupied 3d- 
shell of Co substrate - as in our case, it holds for the Cr 
and O localized magnetic moments. 

As for the difference between the pSIC and the GCA 
results, the pSIC magnetizations are always larger than 
the GCA ones. This is due to the fact that the Hund’s 
rule is strongly obeyed when the electron localization is 
higher like in the pSIG and DFT-I-U case. Moreover, 
for the separate molecules, the magnetizations are local¬ 
ized only at the transition metals and oxygen. In con¬ 
trast, for the periodic systems, the magnetizations are 


also partially localized on N, O, and even carbons of the 
acethylene bridges. 

Even if the magnetization is not localized at the acety¬ 
lene bridges, the covalent bonds between the molecules 
are essential for the spin delocalization within the 
network-building molecules. Without these bonds, the 
molecules would be spin-polarized with the well localized 
integer spin multiplicity, which would not be present at 
the nitrogen atoms. In contrast, for the periodic sys¬ 
tems containing the N atoms and magnetic impurities - 
such as already mentioned diluted magnetic semiconduc¬ 
tor (Ga,Mn)N - the spin polarization is delocalized to the 
nitrogen atoms, with quite substantial contribution!^ 

The more extensive view on the magnetization can be 
derived from the projected densities of states (PDOS) 
onto the 2p-N, 3(i-TM and 2p-0 states, and the total 
DOS, which are presented in Fig. 5; for the results ob¬ 
tained with the pSIG approach. Similar PDOS plots for 
the GGA method are included in Fig. S3 in the supple¬ 
mental information!^ The numbers presented in Table 1 
are obtained from the quadrature of the atomic spin- 
densities presented in Fig. 5. It is in the energy range 
from the bottom of the valence band - not presented in 
the figure - up to the Fermi level. The spin-asymmetry 
of the PDOS is visible for the cases: V, Gr, Mn, Fe, and 
Go without and with the oxygen. The spin-asymmetry of 
the 2p-N PDOS is pronounced even for the cases where 
the Lbwdin spin-polarizations are rather small. Interest¬ 
ingly, for some cases - like Gr - the spin polarization at 
the Fermi level does not show up, but the total contri¬ 
bution to magnetism originates from the deeper states. 
In this case, mostly from Gr and less from the nitrogen 
atoms. 


B. Metallicity 

In Fig. 5, we see that most of P-TM and P-TMO net¬ 
works are metallic. The Fermi levels of some cases - 
namely Cr, MnO, TiO, Ni, Zn - are placed within a little 
energy gap. The half-metallicity is plausible in the cases: 
Ti, Cu, VO, CrO, FeO, CoO, where the more accurate 
GW-calculations would be necessary!^ It is very usual 
for the TM-doped metallic systems, that the Fermi level 
cuts the Sd-states. The effect of such methods like the 
DFT-I-U or the pSIC with respect to the GGA results 
moves the 3d-states away from the Fermi level, if these 
states are not pinned to it. In our systems, however, some 
cases calculated with the GGA approach show clearly the 
TM-based metallicity, while in the pSIC approach, the 
3d states disappear from the Fermi level - e.g. for the 
porphyrin-Mn, -Fe, -VO, -MnO, -CoO and -Ni networks. 
Purely carbon-based metallicities in such magnetic sys¬ 
tems like P-VO, P-Mn and P-FeO networks seem to be 
promissing for spintronic devices, because the magnetic 
couplings in these cases might be long range. Interest¬ 
ingly, there are also the metallic systems where the oxy¬ 
gen states dominate at the Fermi level. These systems 
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FIG. 5. Projected densities of states (PDOS) of the metal- and oxo-metal-porphyrin 2D networks obtained with the pSIC 
approach. Grey color depicts the total DOS, the brown color is the PDOS of 3d-TM states, the blue line is the PDOS of 2p-N 
(summed for 4 atoms), the red line is the PDOS of 2p-0 states. Fermi level is at zero energy. 


are almost half metallic - for instance in the porphyrin- 
NiO, -CuO and -ZnO networks. 


IV. CONCLUSIONS 

2D organic networks are promissing for spintronic 
applications. In this work, we searched for magnets 
among oxo-metal-porphyrins and metal-porphyrins con¬ 
nected with the acetylene bridges. For magnetism, the 
relative energetic positions of the 3d—, 4s— and 4p— 
transition-metal states with respect to the 2p-states of 
the neighbouring atoms are very important. Therefore, 
we used the self-interaction corrected pseudopotential 
scheme (pSIC), within the DFT framework for the cal¬ 
culations of the electronic structure. 

Most of the porphyrin networks with TM and TMO are 
magnetic; except the Sc, TiO, Ni and Zn cases. Addition 
of oxygen increases the magnetizations for the second- 
half 3d-row TMs and the ScO case. The rest of cases: Ti, 
V, Cr and Mn decrease magnetizations after the oxygen 


addition. 

Difference between the total and absolute magnetiza¬ 
tions indicate ferrimagnetism, which is the strongest for 
the porphyrin-V, -VO and -Mn networks. Further exam¬ 
ination of the systems with the Lowdin populations tool 
shows the AF-coupling of the oxygen magnetic moment 
with the TM local moment for the P-VO and P-CrO 
cases, and FM-coupling for the P-FeO, P-CoO, P-NiO, 
P-CuO and P-ZnO systems. Interestingly, in the P-ScO 
case, the magnetic moment is localized at oxygen and 
not at TM, and in the MnO case the situation is oppo¬ 
site. The strongest spin-polarization localized at nitro¬ 
gens, and AF-coupled to the 3d-TM local moment, was 
found for the P-Mn, P-V, P-VO and P-Fe cases. The 
same couplings are ferromagnetic for the P-NiO, P-Cu 
and P-CuO 2D networks. The spin-polarization effects 
are very long range for the porphyrin networks with the 
TMs from V to Co. In these cases, the spin-asymmetry 
of the atomic shells was found very far from the TM 
atoms, and it was well pronounced even on the acetylenic 
bridges. 

In summary, the ferrimagnetic cases found in our 
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studies are characterized by the induced long-range 
spin-polarizations. These systems certainly would be 
good candidates for the high-Tc magnets and will show 
up new properties when are deposited at the surfaces of 
other materials. 
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